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Spatial  heterod3nie  spectroscopy  (SHS)  is  a  Fourier-transform  spectroscopic  technique  that  simulta¬ 
neously  records  all  path  differences  using  a  detector  array.  Compared  to  conventional  Fourier-transform 
spectroscopy  that  measures  interferogram  samples  sequentially  in  the  time  domain,  SHS  is  insensitive 
to  a  changing  scene;  however,  the  effects  caused  by  differences  in  the  detector  elements  and/or  the  optics 
for  each  sample  must  be  addressed  with  a  flatfield  correction.  The  flatfield  correction  is  t3q)ically  a 
characteristic  of  the  instrument  and  does  not  change  with  the  observed  scene.  We  present  three  different 
flatfielding  approaches.  Each  is  based  on  different  assumptions  and  is  applicable  depending  on  the 
instrumental  effects  dominating  the  flatfield.  ©  2006  Optical  Society  of  America 
OCIS  codes:  300.6300,  120.2650,  040.5160. 


1.  Introduction 

Spatial  heterodyne  spectroscopy  (SHS)  is  a  technique 
that  makes  use  of  a  spatially  scanned  interferometer. 
This  class  of  interferometers  is  receiving  increased 
attention  in  the  experimental  spectroscopy  commu¬ 
nity.  ^  The  basic  SHS  configuration  is  similar  to 
that  of  a  Michelson  interferometer  used  for  Fourier- 
transform  spectroscopy  (FTS)  with  the  return  mir¬ 
rors  replaced  by  fixed  diffraction  gratings  as  shown  in 
Fig.  1.  The  simplest  explanation  of  the  SHS  principle 
is  to  consider  the  plane  wavefronts  entering  the  in¬ 
terferometer  as  indicated  by  the  dashed  line  following 
Li  in  Fig.  1.  At  the  exit  of  the  interferometer  the 
resulting  two  wavefronts  are  tilted  with  respect  to 
each  other  with  a  wavenumber-dependent  tilt  angle. 
These  crossed  wavefronts  produce  a  Fizeau  fringe 
pattern  that  is  imaged  on  the  detector  array.  Zero 
spatial  frequency  corresponds  to  the  Littrow  wave- 
number  of  the  gratings.  For  other  wavenumbers  the 
spatial  frequency  is  proportional  to  the  difference  be¬ 
tween  the  signal  wavenumber  and  the  Littrow  wave- 
number.  The  fringe  pattern  recorded  by  the  detector 
array  is  therefore  equivalent  to  the  Fourier  transform 
of  the  incident  spectrum  that  is  heterodyned  around 
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the  Littrow  wavenumber.  A  defining  characteristic  of 
SHS  instruments  is  that  they  achieve  the  theoretical 
resolution  of  the  dispersive  elements  (gratings)  but 
have  the  large  throughput  associated  with  inter¬ 
ferometers  (including  the  possibility  of  field  widen¬ 
ing)  with  no  moving  parts.  A  detailed  description  of 
the  SHS  concept  can  be  found  in  publications  by 
Harlander  et 

One  major  difference  between  SHS  and  a  conven¬ 
tional  FTS  is  that  a  FTS  interferogram  is  recorded  in 
the  time  domain  by  a  single  detector,  whereas  the 
entire  SHS  interferogram  is  recorded  simultaneously 
with  multiple  detector  elements.  As  a  consequence, 
the  SHS  interferogram  is  not  compromised  by  tem¬ 
poral  source  variations,  which  contaminate  the  spec¬ 
tral  information  in  the  FTS  case.  On  the  other  hand, 
every  SHS  interferogram  sample  is  recorded  using  a 
different  detector  element,  and  the  signal  uses  differ¬ 
ent  paths  through  the  interferometer  optics  for  each 
interferogram  sample  (see  Fig.  1).  If  the  sample-to- 
sample  variations  are  not  deemed  negligible,  a  fiat- 
field  correction  must  be  applied.  These  intersample 
variations  are  similar  to  a  temporally  changing  scene 
in  the  FTS  case,  with  the  important  difference  that  in 
SHS  the  variations  are  a  constant  characteristic  of 
the  spectrometer  and  do  not  depend  on  the  detected 
scene,  thus  making  a  scene-independent  correction 
possible. 

We  adapt  the  term  “flatfielding”  for  this  correction 
since  it  is  commonly  used  for  the  elimination  of  pixel- 
to-pixel  variations  in  imaging  detectors,  which  is  a 
closely  related  problem.  Here  we  present  three  SHS 
flatfielding  approaches.  The  first  one  is  suitable  for 
instruments  in  which  the  detector  pixel-to-pixel  sen¬ 
sitivity  variation  dominates.  The  second  one  is  suit- 
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Fig.  1.  Schematic  diagram  of  basic  non-field- widened  SHS  con¬ 
figuration.  The  dashed  lines  illustrate  incoming  wavefronts  and 
the  corresponding  exiting  wavefronts  that  are  crossed  with  an 
angle  of  27.  The  ray  bundle  for  two  interferogram  samples  is  out¬ 
lined  showing  that  only  a  small  section  of  the  interferometer  and 
optics  is  used  for  any  individual  interferogram  sample. 

able  if,  in  addition  to  the  detector  sensitivity  varia¬ 
tions,  the  signal  from  both  interferometer  arms  is  not 
identical,  e.g.,  due  to  the  beam-splitter  performance. 
The  third  approach  is  suitable  for  cases  in  which  a 
spatial  variability  (e.g.,  scratches,  manufacturing  de¬ 
fects,  or  surface  contamination)  in  the  grating  perfor¬ 
mance  also  plays  a  significant  role. 

Before  we  present  the  fiatfielding  techniques,  we 
introduce  the  raw  SHS  interferogram,  which  is  the 
fundamental  equation  for  the  discussion  in  Section  2. 

2.  Raw  Spatial  Heterodyne  Spectroscopy 
Interferogram 

Using  the  small-angle  approximation,  the  ideal  SHS 
interferogram  can  be  written  as^ 

1  r 

I(x)  =  2  (^)[^  cos{27T[4(a  —  ao)tan  0Jx}]da, 

0 

(1) 

where  B^(a)  is  the  wavenumber-dependent  spectral 
density,  uq  is  the  Littrow  wavenumber  for  which  the 
recorded  fringe  frequency  is  zero,  0^  is  the  grating 
angle,  and  x  is  the  pixel  location.  The  leading  factor  of 
V2  accounts  for  the  fact  that,  as  in  FTS,  half  of  the 
incident  signal  leaves  the  interferometer  through  the 
input  port,  and  half  is  detected  at  the  output  port.  To 
simplify  the  following  steps,  we  introduce  the 
heterodyned  fringe  frequency  k  =  4(a  -  ao)tan  0^ 
and  the  heterodyned  spectral  density  B{k)  = 


+  k/4  tan  0J  +  -  k/4  tan  0l)]/4  tan  0^  with  k 

>  0.  Using  these  definitions  we  get 


I(x)  =  2  cos{27TK:r}]dK,  (2) 

0 

showing  that  the  heterodyned  spectral  density  is  es¬ 
sentially  the  Fourier  transform  of  the  modulated  part 
of  the  interferogram.  Note  that  the  signals  from  the 
upper  (a  >  uq)  and  lower  (a  <  uq)  sideband  of  the 
spectral  density  both  result  in  positive  fringe  fre¬ 
quencies  and  therefore  are  superimposed  in  the  spec¬ 
trum  retrieved  by  Fourier  transformation  of  the 
interferogram  unless  one  sideband  is  suppressed  by  a 
spectral  prefilter.  Harlander  et  al.^  detailed  a  tech¬ 
nique  for  unambiguously  separating  these  sidebands 
by  using  a  cross  tilt  in  the  grating  alignment  and 
employing  a  two-dimensional  Fourier  transform  to 
generate  the  spectrum. 

In  practice,  nonideal  instrument  properties  can 
have  additional  effects  on  the  raw  interferogram.  The 
raw  interferogram  measured  along  one  row  of  a  SHS 
detector  array  can  generally  be  represented  as 

I(x)  =  I  S(k){^a^(jc,  k)  +  tB^(x,  k) 

0 

+  2s(x,  K)tj^(x,  K)tB(X,  k) 

X  COSpTTKX  +  @(X,  K)]}dK,  (3) 

where  six,  k)  is  the  modulation  efficiency,  and  @(jc,  k) 
is  a  phase  distortion  term.  The  terms  tA^ix,  k)  and 
ts^ix,  k)  are  the  intensity  transmission  functions  for 
the  signal  passing  through  the  two  interferometer 
arms  A  and  B,  multiplied  by  the  individual  detector 
responsivity.  The  resulting  raw  interferogram  lix)  is 
generally  recorded  in  digital  units  using  an  analog- 
to-digital  converter  (ADU).  The  raw  measurement  is 
subsequently  calibrated  yielding  the  spectral  radi¬ 
ance  of  the  source.  The  calibration  generally  requires 
proper  fiatfielding,  therefore  we  will  use  raw  inter- 
ferograms  in  Section  5. 

Equation  (3)  is  a  general  representation  of  a  SHS 
interferogram  and  allows  for  all  possible  spectral  and 
spatial  dependencies.  For  an  ideal  SHS  instrument, 
tA^(x,  k)  and  ts^ix,  k)  are  identical  (1/4)  and  do  not 
depend  on  x  or  k,  @(jc,  k)  is  zero,  and  six,  k)  is  close  to 
unity  for  zero  optical  path  difference  and  slowly  de¬ 
creases  for  higher  path  differences  due  to  the  finite 
etendue  of  the  interferometer.  ^ 

It  has  been  shown  in  the  work  of  Harlander  et  al.^ 
and  Englert  et  al.^  that  a  phase  distortion  can  be 
corrected  without  impacting  the  envelope  of  the  mod¬ 
ulated  part  of  the  interferogram.  Nevertheless,  we 
will  carry  the  phase  distortion  term  along  to  show 
that  it  does  not  affect  the  fiatfielding  and  that  it  is 
conserved  by  the  fiatfielding.  We  further  make  the 
assumption  that  the  spectral  dependence  of  k) 
and  ts^ix,  k)  is  identical  for  both  arms  so  it  can  be 
separated  into  the  function  i?(K),  and  we  can  write 
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(4a) 


I(x)  = 


B(K)i?(K){fA^(%)  +  tB^(x)  +  2s(x,  k)^a(x)^b(%)COS[2'TT-|OC  +  0(x,  K)]}dK, 


I(X)=  B(K)R(K)[t/(x)  +  tB\x)]dK  + 


2B(K)i?(K)s(X,  K)f^(x)fB(x)C0S[2TTKX  +  ©(x,  K)]dK 


(4b) 


Nonmodulated  term 


Modulated  term 


A  consideration  of  the  effects  contributing  to  the 
spectral  dependences  of  k)  and  ts^ix,  k)  con¬ 
firms  that  this  assumption  is  good:  the  spectral  de¬ 
pendence  due  to  the  spectral  response  of  the  detector 
elements  is  by  design  identical  for  the  signal  of  both 
arms.  The  design  of  both  interferometer  arms  is  iden¬ 
tical  so  that  the  spectral  dependence  of  the  grating 
efficiency,  the  optical  transmittance,  reflectance,  and 
scattering  occurring  on  all  optical  elements  in  the 
instrument  can  safely  be  assumed  to  be  identical  for 
both  arms.  Furthermore,  for  SHS  instruments  that 
have  a  very  narrow  bandpass  [e.g.,  3  nm  in  the  near 
UV  for  the  Spatial  Heterodyne  Imager  for  Meso- 
pheric  Radicals  (SHIMMER)  on  STPSat-1  (Space 
Test  Program  Satellite- 1)^,  spectral  dependencies 
are  not  only  very  similar  for  both  arms,  but  also  very 
small.  Note  that  this  assumption  does  not  imply  a 
50%/50%  beam  splitter,  but  only  that  the  spectral 
dependence  of  the  beam-splitter  reflection  and  the 
transmission  are  equal  for  both  arms. 

3.  Effects  of  Nonideal  Flatfields 

The  goal  of  the  SHS  flatfield  correction  is  to  eliminate 
the  interpixel  variations  introduced  into  the  inter- 
ferogram  through  the  terms  t^ix),  t^ix),  and  s{x,  k). 

The  most  straightforward  nonideal  flatfield  is 
caused  by  a  random  pixel  sensitivity  variation 
across  the  detector,  which  results  in  a  multiplica¬ 
tive  function  applied  to  the  perfect  interferogram 
and  a  decreased  signal-to-noise  ratio  in  the  re¬ 
trieved  spectrum.  If  we  assume  an  otherwise  ideal 
double  sided  interferogram  with  N  pixels  and  a  ran¬ 
dom,  relative  pixel  sensitivity  variation  of  8s/s,  the 
signal-to-noise  ratio  in  the  retrieved  spectrum  will  be 
between  ^N/2s/hs  and  ^j2/Ns/hs,  where  the  two  ex¬ 
treme  cases  result  from  multiplex  noise  propagation 
assuming  a  single,  unresolved  line  and  a  continuum 
spectrum,  respectively.  For  N  =  640  and  a  random 
pixel  sensitivity  variation  of  8s/s  =  5%,  for  example, 
the  uncorrected  flatfield  would  cause  a  signal-to- 
noise  ratio  in  the  spectrum  of  between  357.8  and  1.1, 
depending  on  the  spectral  content  of  the  source.  This 
simple  example  demonstrates  that  flatfield  effects 
can  have  a  serious  impact  on  the  quality  of  the  re¬ 
trieved  spectra,  especially  for  full  spectra  that  do  not 
contain  only  a  few  isolated  emission  lines. 

In  case  the  flatfield  shows  a  significant  structure, 
for  example,  a  scratch  in  the  grating  or  a  significant 


decrease  in  modulation  efficiency  toward  the  edge  of 
the  detector,  it  is  also  appropriate  to  interpret  the 
effect  of  the  flatfield  as  a  change  in  the  instrumental 
line  shape  function.  As  shown  in  Eq.  (4b),  the  uncor¬ 
rected  flatfield  results  in  multiplicative  functions  for 
the  nonmodulated  and  modulated  part  of  the  inter¬ 
ferogram.  The  multiplicative  functions  are  equiva¬ 
lent  to  convolution  kernels  in  the  spectral  domain 
that  have  to  be  applied  to  the  interferogram  terms. 
The  changed  instrumental  line  shape  typically  re¬ 
sults  in  a  lower  spectral  resolution  and  in  systematic 
errors  from  the  incomplete  removal  of  the  contribu¬ 
tions  from  the  nonmodulated  interferogram  term.  In 
general,  flatfield  errors  result  in  higher  noise,  lower 
spectral  resolution,  and  systematic  errors  in  the  re¬ 
trieved  spectrum. 

The  magnitude  of  the  introduced  errors  obviously 
depends  on  the  severity  of  the  flatfield  variations.  In 
Subsection  5.D  we  present  an  example  of  a  measured, 
flatfield  corrected  spectrum  and  the  difference  to 
the  uncorrected  spectrum  by  using  a  monochromatic 
source  (see  Fig.  6). 

4.  Flatfielding  Techniques 

We  consider  three  flatfielding  techniques  in  this  pa¬ 
per,  which  should  be  applied  depending  on  what  in¬ 
strumental  effects  dominate  the  variations  in  the 
flatfield.  In  Subsections  4.A-4.C  we  describe  these 
three  flatfielding  approaches  and  the  conditions  un¬ 
der  which  their  application  is  suitable.  In  Section  5 
we  illustrate  the  flatfielding  problem  using  data 
taken  with  a  breadboard  SHS  instrument. 

A.  Balanced  Arm  Flatfielding  Approach 

The  balanced  arm  flatfielding  approach  is  the  sim¬ 
plest  of  the  three.  It  is  based  on  the  assumption  that 
^a(^)  —  01*  that  the  signal  from  both  interferom¬ 

eter  arms  is  equal  or  balanced.  This  is  a  good  assump¬ 
tion  when  the  spatial  variability  of  and  tsix)  is 
dominated  by  the  sensitivity  of  the  detector  elements. 
In  this  case,  the  flatfield  is  independent  of  the  in¬ 
terferometer  arm  (A  or  B)  that  the  signal  passes 
through.  This  approach  does  not  correct  for  de¬ 
creased  fringe  modulation  due  to  an  efficiency  mis¬ 
match  between  the  interferometer  arms. 

The  idea  behind  this  approach  is  to  measure  the 
nonmodulated  part  of  the  interferogram,  which  can 
be  accomplished  in  a  number  of  ways.  One  way  is  to 
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take  two  measurements,  each  with  one  interferom¬ 
eter  arm  blocked  by  inserting  an  opaque  material 
between  the  beam  splitter  and  the  grating.  The  sum 
of  the  measurement  with  arm  A  blocked  and  the 
measurement  with  arm  B  blocked  yields  the  non- 
modulated  part  of  the  interferogram.  Another 
method  is  to  change  the  optical  path  in  one  of  the 
arms  during  an  exposure  period  by  significantly 
more  than  one  half  of  the  longest  observed  wave¬ 
length,  so  that  the  modulation  is  completely  washed 
out  for  all  wavelengths.  This  can  be  done,  for  exam¬ 
ple,  by  pressure  scanning  one  arm  or  tilting  a  plane- 
parallel  plate  in  one  arm. 

The  nonmodulated  part  of  the  interferogram  can  be 
written  as 


Iu{x)  =  B„(K)i?(K)[f/(x)  +  fB^(x)]dK 


dimensional  transform  to  generate  the  spectrum. ^ 
In  these  cases  the  additive  term  C2  at  zero  spatial 
frequency  can  be  ignored. 

After  phase  distortion  correction,^  6  the  second 
term  of  Eq.  (5)  is  simply  the  cosine  transform  of  the 
spectral  intensity  B{k)  multiplied  by  R{k)/Ci  and 
s{x,  k).  R{k)/Ci  represents  the  spectral  response  of 
the  spectrometer.  According  to  the  convolution  theo¬ 
rem,  the  modulation  efficiency  s(x,  k)  represents  the 
Fourier  transform  of  the  instrumental  line  shape 
function,  which  can  be  verified  by  measuring  line 
sources. 

B.  Unbalanced  Arm  Flatfielding  Approach 

The  unbalanced  arm  flatfielding  approach  does  not 
require  the  balanced  arm  assumption  t^ix)  =  t^ix), 
but  requires  instead  the  less  stringent  assumption 
that  the  integral  in  Eq.  (6)  is  zero  when  averaged  over 
all  interferogram  samples: 


=  \tA{x)  +  tB^{x)\  B„(K)i?(K)dK 

0 

=  C^\tA{x)  +  (5) 

where  Ci  is  a  scalar  that  depends  on  the  spectrum  of 
the  observed  source  5a(x).  The  ratio  of  the  interfero¬ 
gram  in  Eq.  (4b)  and  /[/(x)  yields 


=  C2  +  ^J  B{k)R{k) 

r  tA{x)tB{x) 

X  J  28(X,  k)  7  2,.  ,  ,  2 7" 


tA  (X)  +  ts  (X) 
X  COSpTTKJC  +  @(JC,  k)]  IdK, 


(6) 


where  the  scalar  C2  is  unity  if  the  same  source  is  used 
for  both  measurements.  Using  the  above-mentioned 
assumption  =  tsix)  we  can  write 


(•<)| 


B(K)i?(K)-j  2e(x,  k) 


tA(x)tB(x) 
tA{x)  +  tB^{x) 


X  COSp'TTO:  +  @(JC,  k)] 


=  0. 


(8) 


For  slowly  varying  t{x)  and  e{x,  k)  and  a  negligible 
signal  at  low  spatial  frequencies  this  is  a  good  as¬ 
sumption.  This  approach  allows  the  correction  of  the 
interferogram  if  in  addition  to  pixel-to-pixel  sensitiv¬ 
ity  variations,  the  signal  from  both  arms  is  not  iden¬ 
tical,  e.g.,  because  of  the  beam-splitter  properties. 
Like  the  first  approach,  it  also  does  not  correct  for 
modulation  efficiency. 

The  idea  behind  this  approach  is  to  measure  the 
nonmodulated  part  of  the  interferogram  for  the  sig¬ 
nal  from  each  interferometer  arm  separately.  By 
blocking  interferometer  arm  B  and  A,  one  at  a  time, 
as  described  in  Subsection  4.A  we  get  the  following 
two  interferograms: 


/(jc)  1  f 

Iu(xj  ^  Cl  J  {'S(k)'R(k)s(x,  k) 

X  COSpTTIcr  +  @(x,  K)]}dK.  (7) 

We  have  eliminated  the  interferogram  sample  (jc)  de¬ 
pendence  in  the  nonmodulated  term  and  only  the 
scalar  C2  remains.  After  Fourier  transformation 
into  the  spectral  domain  this  scalar  will  contribute 
solely  to  the  retrieved  spectral  density  at  the  Lit- 
trow  wavenumber  (zero  spatial  frequency).  Since  all 
wavenumbers  contribute  a  signal  at  zero  spatial 
frequency  through  the  nonmodulated  term  in  Eq. 
(4b),  most  SHS  instruments  have  eliminated  any 
real  spectral  signal  at  zero  spatial  frequency  by 
either  choosing  a  spectral  prefilter  that  is  opaque  at 
the  Littrow  wavenumber  or  by  using  a  cross  tilt 
in  the  grating  alignment  and  employing  a  two- 


Ia(x)=  Ba(K)R(K)tA^(x)dK 

0 

=  tA\x)  S,(K)/?(K)dK  =  Ci7(^),  (9) 

0 


7b(x)  =  S„(K)/?(K)^B^(^)dK 

0 

=  tB^{x)\  S,(K)i?(K)dK  =  Ci^B7)-  (10) 

0 

Again,  Ci  is  a  scalar  that  depends  on  the  spectrum  of 
the  observed  source  B^(k).  The  ratio  of  the  interfero¬ 
gram  in  Eq.  (4b)  and  the  sum  of  Ia{x)  and  Ib{x)  yields 
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,00 


(11) 


Ia(x)  +  hix) 


1 


tA(x)tB(x) 

B(^k)R(^k)2s(X,  ^  ^C0S[27TKJC  +  @(JC,  K)]dK. 


Modulated  term 


The  next  step  is  to  subtract  the  scalar  C2.  If  the  same 
source  is  used  for  the  flatfield  measurements  Ia(x) 
and  Ib{x)  and  the  interferogram  I{x),  then  C2  is  unity 
and  we  can  eliminate  it  by  subtracting  unity  from 
Eq.  (11).  In  the  general  case  in  which  a  different 
source  is  used  for  the  interferogram  I(x),  we  can  use 
the  assumption  of  Eq.  (8)  and  determine  C2  by  calcu¬ 
lating  the  average  of  the  right-hand  side  of  Eq.  (11). 
After  the  subtraction  of  C2  we  can  divide  the  modu¬ 
lated  term  in  Eq.  (11)  by  {2[Ia{x)Ib{x)Y^^}/[Ia{x)  + 
Ib(x)]  to  get  the  flatfield  corrected  interferogram  Ic(x): 


lc{x)  = 


5(K)i?(K)8(jC,  k)COS[27TKJC  +  @(JC,  K)]dK. 


(12) 


This  result  is  similar  to  the  balanced-arm  approach, 
but  it  properly  corrects  for  the  effect  of  unequal  or 
unbalanced  signal  levels  from  the  interferometer 
arms. 


C.  Phase-Shift  Flatfielding  Approach 

The  phase-shift  flatfielding  approach  is  motivated 
by  measurements  with  spectrometers  that  show 
grating  imperfections  such  as  manufacturing  de¬ 
fects,  scratches,  or  surface  dust  particle  contamina¬ 
tion.  These  defects  generally  increase  the  scattering 
at  the  affected  grating  location  and/or  decrease  the 
grating  efficiency.  The  affected  interferogram  sam¬ 
ples  typically  show  lower  overall  signal  levels  and 
decreased  modulation  efficiency.  These  effects  are 
generally  wavelength  independent,  especially  for 
narrowband  instruments. 

The  goal  of  this  flatfielding  approach  is  to  quantify 
the  modulated  and  nonmodulated  contribution  to  a 
monochromatic  interferogram  to  provide  a  general 
flatfield  correction  that  includes  the  effect  of  grating 
imperfections.  This  approach  corrects  for  all  flatfield 
effects,  but  it  is  also  the  most  complex.  It  requires  a 
monochromatic  light  source  and  the  ability  to  change 
the  optical  path  length  in  one  interferometer  arm. 

The  approach  requires  the  measurement  of  three 
interferograms,  each  using  the  same  monochromatic 
source,  but  different  optical  path  lengths  in  one  of  the 
interferometer  arms.  The  change  in  optical  path 
length  can  be  achieved  by  moving  one  grating,  pres¬ 
sure  scanning  one  arm,  or  tilting  a  plane-parallel 
plate  in  one  arm.  The  resulting  interferograms  can  be 
written  as 


Ii{x)  =  S(Ko)-R(Ko){^A^(a:)  +  tB^{x)  +  2e{x,  Ko)^a(^) 

X  fB(%)COS[2TTKoX  +  @(X,  k)  +  (p;]} 

=  N(X,  Ko)  M(x,  Ko)COS[27TKoX  +  @(x,  k)  +  cpj], 

i  =  1,  2,  3  A  cpi  cp2  ^  93-  (13) 

The  next  step  is  to  determine  the  phases  Oj(jc)  = 
[27tkoX  +  0(jc,  k)  +  cpj  for  each  measurement  (i  = 
1,  2,  3)  and  each  interferogram  sample  x.  This  can  be 
achieved  by  using  the  technique  described  by 
Harlander  et  al.^  and  Englert  et  al.^  It  consists  of 
Fourier  transformation  of  the  interferogram,  isola¬ 
tion  of  the  spectral  line  for  positive  frequencies,  and 
backtransformation  into  the  interferogram  domain, 
where  the  ratio  of  the  imaginary  part  and  the  real 
part  is  equivalent  to  the  tangent  of  the  phase.  Note 
that  this  technique  to  determine  the  phase  also  ac¬ 
counts  for  phase  distortions  and  slight  changes  in  kq 
that  might  result  from  a  grating  tilt  during  the  trans¬ 
lation  of  a  grating.  After  determining  the  phases 
fi>j(jc),  we  can  calculate  the  nonmodulated  contribu¬ 
tion  N{x,  Ko)  and  the  amplitude  of  the  modulated  con¬ 
tribution  M{x,  Ko)  of  the  interferograms  for  each 
interferogram  sample  using  the  above  measure¬ 
ments: 


I Ax)  - 1 Ax) 

^  cos  ^^(x)  -  cos  ^ 

N{x,  kq)=Ii{x)—M{x,  ko)cos  fi>i(jc). 


(14) 

(15) 


where  i  and  j  in  Eq.  (14)  are  chosen  for  each  inter¬ 
ferogram  sample  so  that  the  denominator  is  never 
zero.  This  is  the  reason  why  three  measurements 
with  different  phase  offsets  are  necessary  and  two  are 
not  sufficient.  It  is  best  to  choose  i  and  j  so  that  the 
absolute  value  of  the  denominator  is  as  large  as  pos¬ 
sible  to  avoid  the  division  by  a  small  number,  which 
increases  the  uncertainty  of  the  result  when  noise  is 
present.  A  good  choice  for  the  three  phase  offsets 
9i/2/3  is,  for  example,  0°,  90°,  and  180°,  which  causes 
the  absolute  value  of  the  denominator  in  Eq.  (14)  to 
fall  always  between  approximately  1  and  yj2. 

We  can  now  eliminate  the  nonmodulated  contribu¬ 
tion  of  any  interferogram  by  dividing  the  interfero¬ 
gram  by  ^(jc,  Ko)  and  subtracting  the  mean,  similar  to 
the  unbalanced  arm  approach.  Dividing  the  result  by 
M{x,  Ko)  will  remove  the  variations  in  modulation  ef¬ 
ficiency  and  the  effects  of  the  unbalanced  arm  signals’ 
assumption  of  their  wavelength  independence. 
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Table  1.  Key  Design  Parameters  of  the  Visible/Near-Infrared  SHS 
Breadboard  Instrument 


Parameter 

V  alue/Description 

Resolving  power 

~13  000 

Bandpass 

-700-900  nm 

Beam  splitter 

Cubic,  BK7 

Gratings 

60  grooves/mm,  blazed 

Exit  optics 

Telecentric  lens. 

Magnification:  0.5 

Spectral  filter 

Hoya  R-72 

Detector  array 

Micro-Pix  M640 

640  X  480,  7.4  ijim  pitch 

Light  source 

Digikey  M7805I 

780  nm  laser  diode  module 

5.  Flatfielding  Examples  Using  Measured  Data 

A.  Visible/Near-Infrared  Breadboard  Instrument 

The  data  presented  in  the  following  subsections  were 
measured  with  a  visible/near-infrared  breadboard 
SHS  instrument.^  The  key  instrument  parameters 
and  the  light  source  used  for  the  measurements  are 
summarized  in  Table  1.  Figure  2  shows  the  inter¬ 
ferometer  of  the  breadboard  instrument.  All  the  com¬ 
ponents  of  this  instrument  are  commercially  available. 
The  light  from  the  source  was  diffused  at  a  small  di¬ 
ameter  iris  and  collimated  by  a  spherical  lens  before 
entering  the  interferometer.  The  gratings  are  imaged 
on  the  detector  array.  All  the  measurements  have 
been  dark-field  corrected.  The  two-dimensional  CCD 
image  was  binned  to  form  ten  rows  (interferograms) 


Fig.  2.  Interferometer  of  the  visible/near-infrared  breadboard  in¬ 
strument.  The  input  is  on  the  lower  left-hand  side  of  the  cubic 
beam  splitter.  The  gratings  are  mounted  using  commercial  three- 
point  mounts.  The  grating  on  the  right  is  mounted  on  a  translation 
stage  and  can  be  moved  away  or  toward  the  beam  splitter.  All  the 
components  of  this  breadboard  spectrometer  are  commercially 
available  (see  also  Table  1). 


Fig.  3.  Top  panel:  raw  monochromatic  interferogram  (black)  and 
the  nonmodulated  signal  (gray)  measured  by  blocking  one  inter¬ 
ferometer  arm  at  a  time  and  adding  the  two  measurements.  Center 
panel:  interferogram  after  application  of  the  balanced  arm  flat- 
fielding  method.  Bottom  panel:  real  (black)  and  imaginary  (gray) 
parts  of  the  spectrum  calculated  from  the  corrected  interferogram 
by  phase  correction,  zero  filling,  and  Fourier  transformation.  Since 
an  unresolved  line  source  was  used  for  this  measurement,  the 
spectrum  represents  the  instrumental  line  shape  function. 


with  640  samples  each.  For  the  following  discussion, 
we  present  the  data  from  row  7. 

B.  Example  for  Balanced-Arm  Flatfielding 

The  upper  panel  of  Fig.  3  shows  a  raw,  monochro¬ 
matic  interferogram  (black)  and  the  measured  non¬ 
modulated  component  (gray).  The  measurement  of 
the  nonmodulated  component  was  obtained  by  block¬ 
ing  one  interferometer  arm  at  a  time  by  using  a  piece 
of  black  paper  and  adding  the  two  measurements. 
The  nonmodulated  part  of  the  interferogram  shows 
little  pixel-to-pixel  variation,  indicating  very  similar 
pixel  sensitivity.  A  slight  decrease  in  the  nonmodu¬ 
lated  signal  can  be  seen  at  approximately  pixel  560, 
which  is  also  visible  as  a  signal  decrease  in  the  raw 
interferogram.  The  center  panel  of  Fig.  3  shows  the 
interferogram  after  the  balanced  arm  flatfielding  ap¬ 
proach  was  applied.  Since  the  pixel-to-pixel  sensitiv¬ 
ity  variations  are  small  in  this  example,  the  corrected 
interferogram  looks  similar  to  the  raw  interferogram. 
The  fiatfield  correction  improves  the  problem  at  ap¬ 
proximately  sample  560,  but  even  after  the  correction 
a  slight  irregularity  can  be  identified,  which  indicates 
that  it  is  not  caused  by  a  difference  in  detector  sen¬ 
sitivity.  The  bottom  panel  of  Fig.  3  shows  the  real  and 
imaginary  parts  of  the  spectrum  that  were  retrieved 
by  zero  filling  and  Fourier  transforming  the  interfero¬ 
gram.  The  remaining  instrumental  effect  at  approx¬ 
imately  sample  560  is  so  small  that  the  real  part  of 
the  spectrum  shows  no  significant  deviation  from  a 
sine  function,  which  is  expected  from  the  finite  length 
interferogram  and  an  unresolved  laser  line  source. 
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Fig.  4.  Top  panel:  nonmodulated  interferogram  components  mea¬ 
sured  with  one  arm  blocked.  Center  panel:  correction  factor  for  the 
multiplicative  term  of  the  modulated  interferometer  part  that  is 
caused  by  the  unbalanced  arm  contributions.  Bottom  panel:  cor¬ 
rected  interferogram. 


C.  Example  for  Unbalanced-Arm  Flatfielding 

The  upper  panel  of  Fig.  4  shows  the  nonmodulated 
components  of  the  spectrum  for  the  two  interferom¬ 
eter  arms.  They  are  measured  by  blocking  one  arm  at 
a  time.  The  contributions  from  both  arms  are  not 
perfectly  balanced  for  this  interferometer.  Some 
features  are  common  to  both  arms,  others  are  not. 
For  example,  the  feature  at  pixel  560  is  present  only 
in  one  interferometer  arm,  indicating  that  it  is  a 
scratch,  a  manufacturing  defect,  or  contamination  on 
the  corresponding  grating. 

The  middle  panel  of  Fig.  4  shows  the  term 
{2[/a(^)^5(^)]^^^}/[^a(^)  +  hi^)]  that  is  needed  to  cor¬ 
rect  for  the  effect  of  the  unbalanced-arm  contribu¬ 
tions  on  the  envelope  of  the  modulated  part  of  the 
interferogram  [Eqs.  (11)  and  (12)].  Except  for  the  re¬ 
gion  at  approximately  pixel  560,  this  correction  factor 
is  always  within  0.1%  of  unity  for  this  example. 

The  lower  panel  of  Fig.  4  shows  the  interferogram 
after  application  of  the  unbalanced  arm  flatfielding 
method. 

D.  Example  for  Phase-Shift  Flatfielding 

Figure  5  shows  the  measurements  Ii(x),  l2(x),  and 
Is(x)  [Eq.  (13)]  around  the  area  of  the  grating  imper¬ 
fection.  The  different  phase  shifts  of  the  interfero¬ 
gram  were  achieved  by  moving  one  grating  by  using 
a  micrometer  driven  translation  stage.  The  middle 
panel  shows  the  modulated  and  nonmodulated  con¬ 
tributions  to  the  interferogram  as  calculated  using 
Eqs.  (14)  and  (15).  The  lower  panel  shows  the  cor¬ 
rected  interferogram.  The  improvement  made  by  this 
flatfielding  approach  to  the  area  at  approximately 
pixel  560  can  be  seen  clearly  by  comparing  the  raw 
interferograms  in  the  upper  panel  to  the  corrected 
interferogram  in  the  lower  panel. 

Figure  6  shows  the  nonmodulated  and  modulated 
components  for  the  entire  interferogram  in  the  upper 
panel  and  the  corrected  interferogram  in  the  center. 


Pixel  Number 


Fig.  5.  Top  panel:  raw  interferograms  for  three  different  grating 
locations  resulting  in  a  phase  shift  for  each  interferogram.  Center 
panel:  nonmodulated  interferogram  part  (black)  and  the  envelope 
of  the  modulated  interferogram  part  (gray)  as  calculated  from  the 
raw  interferograms.  Bottom  panel:  corrected  interferogram. 


The  lower  panel  shows  the  retrieved  spectrum  and 
the  difference  between  the  corrected  and  the  uncor¬ 
rected  spectrum  multiplied  by  10. 

The  example  data  we  used  to  illustrate  the  flat¬ 
fielding  techniques  are  of  high  quality,  i.e.,  the  beam 
splitter,  grating,  exit  optics,  and  detector  elements 
did  not  cause  major  flatfield  inhomogeneities.  De¬ 
pending  on  the  wavelength  region,  such  high  quality 
components  might  not  always  be  available.  In  these 
cases,  a  proper  flatfield  correction  gains  in  impor¬ 
tance. 


Fig.  6.  Top  panel:  nonmodulated  interferogram  part  (black)  and 
the  envelope  of  the  modulated  interferogram  part  (gray)  as  in  the 
center  panel  of  Fig.  5  but  for  the  entire  interferogram.  Center 
panel:  corrected  interferogram.  Bottom  panel:  real  part  of  the  cor¬ 
rected  spectrum  and  the  difference  between  the  corrected  and 
uncorrected  spectrum  multiplied  by  10. 
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6.  Conclusion 

Variations  in  optics  and  detector  pixels  have  long 
been  recognized  as  a  data  analysis  issue  for  all  spec¬ 
trometers  that  use  detector  arrays.  We  have  pre¬ 
sented  three  methods  to  correct  for  pixel-to-pixel 
variations  in  spatial  heterodyne  spectroscopy.  When 
the  flatfield  error  is  dominated  by  the  sensitivity  vari¬ 
ations  in  the  detector  array  the  balanced  arm  ap¬ 
proach  is  appropriate.  When  significant  differences  in 
the  signal  from  the  two  interferometer  arms  exist,  the 
unbalanced  arm  approach  is  appropriate.  When  the 
gratings  exhibit  significant  irregularities  such  as 
scratches,  the  phase-shift  fiatfielding  approach  is  ap¬ 
propriate.  All  fiatfielding  corrections  require  addi¬ 
tional  measurements.  The  easiest  way  to  obtain  the 
necessary  fiatfield  measurements  for  the  first  two 
methods  is  to  block  one  interferometer  arm  at  a  time. 
The  third  method  requires  a  change  of  the  optical 
path  length  in  one  arm,  which  is  usually  harder  to 
achieve,  but  is  possible  by  moving  the  grating,  pres¬ 
sure  scanning  one  arm,  or  tilting  a  plane-parallel 
plate  in  one  arm. 

In  principle,  any  interferogram  measured  by  one 
particular  spectrometer  can  be  fiatfield  corrected  us¬ 
ing  one  set  of  fiatfield  measurements.  However,  care 
needs  to  be  taken  to  eliminate  effects  that  could 
change  the  fiatfield  such  as  thermal  distortion  of  the 
instrument  causing  a  shift  or  magnification  change  of 
the  imaged  interferogram.  If  the  fiatfield  is  likely  to 
change  with  time,  e.g.,  due  to  extreme  thermal  con¬ 
ditions,  it  may  be  necessary  to  record  fiatfield  mea¬ 
surements  periodically  during  the  observation. 
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